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Of all the fundamental parameters of the Sun (diameter, mass, temperature...), the grav- 
itational multipole moments (of degree I and order m) that determine the solar moments of 
ff^ ■ inertia, are still poorly known. However, at the first order (1=2), the quadrupole moment is 

^ ' relevant to many astrophysical applications. It indeed contributes to the relativistic perihelion 

advance of planets, together with the post-Newtonian (PN) parameters ; or to the precession 
' of the orbital plane about the Sun polar axis, the latter being unaffected by the purely rela- 

. tivistic PN contribution. Hence, a precise knowledge of the quadrupole moment is necessary 

' for accurate orbit determination, and alternatively, to obtain constraints on the PN param- 

■ eters. Moreover, the successive gravitational multipole moments have a physical meaning : 

they describe deviations from a purely spherical mass distribution. Thus, their precise deter- 
mination gives indications on the solar internal structure. Here, we explain why it is difficult 
Oh' to compute these parameters, how to derive the best values, and how they will be determined 

^ . in a near future by means of space experiments. 

t-^ ' Keywords : sun ; solar rotation ; solar gravitational moments ; 

. general relativity ; perihelion advances of planets. 
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■ 1 Introduction 



The study of the rotation of stars is not trivial. In theory, the problem is 
exceedingly simple and can be formulated as follows. Consider a single star that 
rotates about a fixed direction in space, with an angular velocity fi. Let us first 
assume that, for 17 = 0, the star is a gaseous body in gravitational equilibrium. 
The problem is to determine the outer shape of the star when the initial sphere 
is set rotating at an angular velocity ft. Such studies were conducted for the 
first time by Milne (1923), then fully achieved by Chandrasekhar (1933). 

The second point which arises is to ask oneself what will happen if 17 is 
not constant, not only in latitude (differential rotation) but also throughout 
the body, from the surface to the core. Today, astronomers are faced with such 
problems, not only in the solar case, but also for stars. With the advent of 
sophisticated techniques such as interferometry, one is now able to accurately 
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determine the geometrical shape of the free boundary of stars, such as Altair or 
Achernar for which observations of the geometrical enveloppe have been made 
by Belle et al. (2001), and Domiciano de Souza (2002). But it would be of little 
or no interest to observe the geometric shape of a star if one would not be able to 
infer some information on stehar physics. With such an approach, the purpose of 
theoreticians is enumerate all the possible angular velocity distributions (from 
the center to the surface) that are compatible with the observed stellar surface. 
For stars, Maeder (1999) examined the effects of rotation and wrote the equation 
of the surface with a rotation law which is differential, but only in the surface 
layer. 

In other words, the knowledge of the angular velocity distribution from the 
core to the surface, together with the knowledge of the density function (related 
to the pressure function), completely determine the outer shape of the stars. Dif- 
ferent techniques exist to observe such a figure. Once accurately determined, one 
would be able to go back to the physical properties of the body. This approach 
is called "Theory of Figures". 

The Theory of Figures has been widely used in geophysics and is still used in 
specific cases, such as for the planet Mars, with an incredible accuracy ( J2 Mars 
= 1.860718 10"'^ according to Yuan et al. (2002), from a 75th degree and or- 
der model). Since the pioneering work of Clairaut (1743), who was the first to 
compute the flattening of a rotating body, and Bruns (1878) who introduced 
the concept of "Figure of the Earth", considerable work has been done, mainly 
by Radau (1885), Wavre (1932), Molodensky (1988), Moritz (1990), and many 
others since then. The basic principle of this theory is to determine the outer 
surface of a rotating body, assuming a knowledge of the mass and angular ve- 
locity distributions. In the case of the Earth, the main emphasis was on the 
determination of a global ellipsoid. The theory of hydrostatic equilibrium was 
considered as the best means for determining the terrestrial flattening /, un- 
til the advent of artificial satellites in 1957. Today, considerable accuracy is 
reached that allows the inverse problem to be solved, that is to say to determine 
anomalies inside the Earth mantle. The Earth external gravitational potential 
is developed in a serie of so-called spherical harmonics of degree I and order 
m. The dimensionless coefficients represent the different distortions from a pure 
sphere. Thanks to artificial satellites, they are now computed up to a very high 
degree (J2 Earth = 4.84165209 lO""* according to the EGM96 tide free model 
adopted by the International Earth Rotation Service in (McCarthy, 2003) , from 
a 360th degree and order model), and substantial progresses are made in the 
knowledge of our planet. 

Curiously, for the Sun, the situation is hardly that of the Earth which pre- 
vailed before 1957. The flattening of the Sun is stih poorly known and we hardly 
know if there is a variability hnked to the solar cycle. The Theory of Figures 
was not developed for the Sun until Rozelot et al. (2001). The difference be- 
tween the Sun and the Earth lies in the fact that, in the latter case, the rotation 
can be considered as solid and the outer shape is mainly due to strong density 
variations inside the body, on which are superimposed tides, oceanic currents 
and atmospheric effects. On the contrary, in the solar case, the density can be 
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Fig. 1 - Left : 3d sketch of the {1= 2, m= 0) mode; the Sun is nearly a spheroid. 
Right : 3d sketch of the (/= 10, m= 6) mode ; distortions from a spherical surface 
appear, on a very exagerated scale. 

considered as nearly homogeneous on successive shells, whereas the rotation is 
strongly dependent on the latitude. Solar modeling is further compHcated by 
the fact that this differential rotation is not only a surface phenomena, as can 
be seen for example by observing spots or faculae, but is anchored in depth, 
at least down to 0.7 solar radii. The computation of the solar gravitational po- 
tential may be not too difficult in its formalism, and is nearly the same as for 
the Earth itself. However the computation of the centrifugal potential is very 
complicated and it cannot be reduced to computation of potentials on successive 
cylinders (or thin zonal rings) in which the rotation is taken as uniform. More- 
over, we have shown that the centrifugal force corresponding to the commonly 
adopted differential rotation law (a development in cosine of the colatitude) 
does not derive from a potential. Hence, we proposed an alternative rotation 
law (Rozelot and Lefebvre, 2003). The complexity of the rotation profile (see 
for example Eq. indicates that the photospheric shape is highly sensitive 
to the interior structure. Thus, in principle, accurate measurements of the limb 
shape distortions, which we call "asphericities" (i.e. departures from the "he- 
Hoid", the reference equilibrium surface of the Sun) , combined with an accurate 
determination of the solar rotation provides useful constraints on the internal 
layers of the Sun (density, shear zones, surface circulation of the plasma...). Fig. 
n] shows such asphericities that can be seen at a given spatial resolution. Al- 
ternatively, theoretical upperbounds could be inferred for the flattening which 
may exclude incorrect /biased observations. 

Thus, even if we know that the Solar Theory of Figures is limited due to trun- 
cation errors, the learning is rich. 

In this paper, we will successively focus on : (i) the definition of the solar 
gravitational moments, (ii) how to compute these moments, (iii) their relevance 
to astrophysics and (iv) how to measure such quantities. To conclude, we will 
show why space missions are needed. 
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2 Definition of the solar gravitational multipole 
moments 



For an axially symmetric distribution of rotating matter, the outer gravita- 
tional field can be expressed as : 



r 



(1) 



where G is the gravitational constant ; Mq and Rq , the solar mass and radius ; 
J2n are the gravitational multiple moments; P2n, the Legendre polynomials; 
r and 9, respectively the distance from the Sun centre and the angle to the 
symmetry axis (colatitude). 
The first terms are called 

- for n = 1, J2 the solar quadrupole moment 

- for n = 2, J4 the solar octopole moment 

- for n = 3, Je the solar dodecapole moment 

For the Sun, as we will see later, terms of higher degree are nearly null. The J„ 
coefficients are dimensionless quantities providing information on how the mass 
and velocity distributions act inside the Sun to finally render the outer visible 
shape non spherical. The amplitude, peak to peak, of such asphericities does not 
exceed some 20 milliarcseconds (abbreviated as "mas"), as shown by the theory 
itself (Lefebvre and Rozelot, 2004) or by observations (Rozelot et al, 2003). 



3 Computation of the solar gravitational moments 

The starting point is to consider the equation of motion for an ideal fiuid 

px = pg- Vp (2) 

where x is the acceleration of the fluid ; p is the density ; g, the gravitational 
force ; and Vp, the pressure gradient. For hydrostatic equilibrium, there is no 
motion, hence x = 0. For a nonrotating body, the force acting on particles is 
given by 5 = —Vcfigrav, where ipgrav is the gravitational potential. Hence Eq. [3 
reduces to ^ 

= -pV(l)grav - Vp (3) 

When the body is rotating, the force acting on particles is the gradient of the 
gravity potential {iptotai = 4>grav+4>rotation) ; that is, we must add the centrifugal 

^The curl of this equation is Vp X ^4> = Oi so that the normals to surfaces of constant p 
and 4) point in the same direction. As a consequence, the surfaces of constant p and 4> coincide. 
Using the perfect gas equation of state (with constant chemical composition), it can be shown 
that surfaces of constant p, (p, T and P all coincide. This is known as the Von Zeipel theorem : 
any internal source of distortion in the gravitational field at the surface will manifest itself as 
a change of shape in the solar surface layer. Thus, measuring the shape of the surface layers 
is equivalent to measuring surfaces of constant gravitational potential. 
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force, so that Eq. El becomes 

Vp = -pV(t)grav + 9) s (4) 

where fl{r, 9) is the angular velocity and s a vector perpendicular to the rotation 
axis directed outwards. 

There are two ways to compute the coefficients J„. 

• The first method solves Eq. by considering that asphericities are small 
quantities : the solution is the sum of a spherically symmetric potential (indexed 
by a subscript 0) and a perturbation (nonspherically symmetric, indexed by a 
subscript 1) : 

4>grav{r, 9) = (j)grav o(f ) + (f)grav l{r^ 9) (5) 

Such an approach was taken by Goldreich and Schubert (1968), Ulrich and 
Hawkings (1981), Pijpers (1998), Armstrong and Kuhn (1999), Godier and 
Rozelot (1999) or Roxburgh (2001). Differences lie in the way Eq. 01 is pro- 
jected on appropriate coordinates. Except for Armstrong and Kuhn who used 
a series of vector spherical harmonics to avoid truncation errors, all the above 
mentioned authors made a projection on Legendre polynomials P2n{cos9), so 
that 

oc 

grav l,2n {r) P2n{cos9) 

n=l 

According to this, Eq. |2 supplemented with Poisson's equation 

'^'^(l^grav = 'inGp (6) 

and boundary conditions at r = 0, where (pgrav = or at the surface r= Rq 
where (pgrav = ^grow out , finally lead to 

J2n — Qj^^ ^ (t>grav l,2n(^o) (7) 

The function (jigrav i,2n is the solution to a differential equation requiring the 
knowledge of p(r) and fl{r, 9). Different models can be used for the density or the 
rotation law. In general p{r) is taken either from Richard (1999), Christensen- 
Dalsgaard et al. (1996) or Morel (1997). The rotation law can be either the 
standard rotation law (in cosine, see Eq.jHjl or can be derived from the analyti- 
cal law of Kosovichev (1998) or Corbard et al. (2002) (see Eq. ??). Fig.El-left- 
shows solar rotation profiles, from the surface down to 0.3 Rq, obtained by a 
1.5 dimensional inversion of data from the MDI (Michelson Doppler Imager) 
experiment on board the SOHO mission (Di Mauro, 2003). Fig.[2l-right- shows 
the Corbard rotation model deduced from observations of the MDI /-modes 
between May 1996 and April 2001 (Dikpati et al, 2002). 

A complete expression oi4>grav 1,2 and 4>grav 1,4 was provided by Armstrong and 
Kuhn (1999), using the standard rotation law. However, surface plasma observa- 
tions allow to constrain analytical rotation models. The first attempt to derive 
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an analytical rotation law from helioseismic data has been made by Kosovichev 
(1996). Using this law, Godier and Rozelot (1999), then Roxburgh (2001) com- 
puted the gravitational moments (see Table QJ. The discrepancy between the 
values obtained can be explained by the use of different density models (the 
model of density used here is that of Richard, 1999) and by the way the differ- 
ential equation Eq. Ejl is integrated. Kosovichev noted that a subsurface shear 
layer results when the helioseismically obtained internal rotation is matched 
with the surface rotation. Hence, we suspect that the subsurface rotation rate, 
in this very thin layer (called the leptocUne), may play an important role. To 
check this point, we used the rotation model described in (Dikpati et al, 2002). 
This model presents the following characteristics : (i) the rotation rate is con- 
stant in the radiative interior (core), i.e. no{r,0) — 435 nHz as in Kosovichev 
(1996) ; (ii) the location and width of the tachocline are assumed to be inde- 
pendent from latitude ; (iii) the rotation rate at the top of the tachocHne and 
at the surface are given by 

n{rcz,d) = ^eq{r) + a2Cos^9 + a^cos'^O (8) 

with 02 = -61 nHz and 04 = -73.5 nHz ; (iv) the radial rotation gradient 
near the surface is assumed to be constant at a given latitude and its latitu- 
dinal dependence is given by 13{6), which is a polynomial function of cosO : 
P{6) = Pq + (3y,{cos9Y + f3e{cos9)^. The transition between layers (of constant 
dfl/dr at a given latitude) is constructed using error functions, i.e. Sa;(r) = 
0.5 {1 -|- er/[2(r — rx)/'Wx]} where w is the width of the transition, and x stands 
either for the tachocHne (tac), the convective zone (cz) or the surface (s). The 
final expression of the solar internal rotation law is given by Eqs. 1, 2 and 3 of 
Dikpati et al, (2002) or, alternatively, Eqs. 1, 2, 4 and 5 of (Corbard et al, 2002). 
Values of the parameters flo, 02 and 04 lead to flcz = 453.5 nHz (but other 
values of 02 and 04 are found by other authors -see Lefebvre and Rozelot, 2004) ; 
we adopted wtac = Wcz = Ws = 0.05 Rq, rtac = 0.69 Rq, rcz = 0.71 Rq and 
Ts = 0.97 Rq. Parameter /3o has been taken as 437 nHz/RQ, instead of 891.5 
used by Kosovichev(1996), and we used = -214 and = -503 {nHz/RQ^. 
The main advantage of this formalism resides in a modeling of the radial dif- 
ferential rotation in the subsurface. The order of magnitude of this outward 
gradient is —5.7 10"^^ m~^s~^ at the equator. In a shell with a thickness of 
about 0.03 Rq, the radial rotation gradient dfl/dr is negative from the equator 
to about 50° of latitude, crosses zero around 50° and is positive beyond this 
latitude. Basu et al (1999) deduced a similar behaviour from their analysis of 
the splitting of high-degree /-modes, but finding a reversal of the radial gradient 
in a zone above 0.994 Rq only. 

After solving equation Eq. it turns out that 

J2 = 2.28 10-^ 

This computed value is very sensitive to the surface parameters (02 and 04, see 
Eq.lHl (latitudinal rotation gradient), and to the latitudinal differential rotation 

^Dikpati et al. (2002) alternatively considered /3-j = and ft = -1445 (uHz/Rq), but 
computations do not show strong differences in J2, only 3 per 1000. 
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Fig. 2 — Left : Solar rotation profiles from 0.3 Rq to the surface for different latitudes, 
from the Equator (top curve) to the Pole (bottom curve), obtained by a 1.5 dimensional 
inversion of SOHO-MDJ data. Courtesy of Maria di Mauro, 2003. Center : Radial 
gradient of the angular velocity as a function of latitude. See (Corbard et al., 2002) for 
further details. Right : Solar rotation profiles from 0.55 Rq to the surface for different 
latitudes, computed every 10° from the Equator (top curve) to the Pole (bottom 
curve), with the analytical rotation model of Corbard et al. (2002). This model takes 
into account the sub-surface gradient, as shown in the center plot (dotted line). 



in the convective zone. It can reach a maximum of 2.62 10~^ and a minimum 
of 1.80 10^'', that is to say a 15% variation. The situation is worse for J4 where 
the excursion of the computed values is more than 25%. We conclude that (i) 
helioseismic rotation rates lead to underestimate values of J„ by comparison to 
values deduced from the theory of Figures (see next section and Table and 
(ii) the multipole moments (J2 and J4), seem to be sensitive to the physical 
mechanisms which act at the near surface. Mecheri et al. (2004), using the same 
formalism but a different model of density, reached nearly the same conclusions. 

• The second method consists in writing the radius r as a function of the colat- 
itude 6, 



r = Rsp 



1 + ^ C2nP2nicOs6 



(9) 



where Rgp is the radius of the best fitting sphere passing through the poles and 
the equator of the Sun (for an ellipsoid of revolution with an equatorial radius a 
and a polar radius &, Rsp — (a^6)^/'^). It is only a matter of algebra to substitute 
the above development (Eq. |2l in the external gravitational potential (Eq. . It 
is a bit more complicated to determine a centrifugal force in Eq.01corresponding 
to the rotation law in (r, 9) and which derives from a potential. Adopting the 
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following rotation law fulfills this requirement : 



oo 



-I 1/2 



(10) 



where Vp = r/Rsp is the normaHzed radial distance from the solar center at a 
given depth and Clpoi is the angular velocity at the pole. Using the solar Green- 
wich data base, Javaraiah and Rozelot (2002) fitted the parameters : 

62 = +0.442 and 64 = +0.056 at the surface {vp = 1) ; flpoi = 381.8 nHz 

The resulting centrifugal potential can then be expressed on the basis of Legen- 

dre polynomials and the total potential is written as 4>totai = ^poi R% h{f, Jn) 
where / is the flattening of the Sun {f = {a — b)/a with the above definitions). 
To determine the J„, it is sufficient to consider (ptotai as independent from 6, 
that is, the coefficients of P2n vanish. The fact that the analytical differential 
rotation models lead to a significantly lower value of the flatness in compari- 
son to a uniform rotation model can only be interpreted in terms of a positive 
outward rotation gradient in the subsurface. 

A full description of the method can be found in Rozelot and Lefebvre (2003) 
or Lefebvre (2003). 

The two methods are slightly different. The flrst one provides the coefficients 
Jn if the mass distribution and rotation rate are known. It does not yield a de- 
termination of the flattening /, and there is no means to link / with J„ without 
strong assumptions. By contrast, the second method requires an analytical ro- 
tation law (deriving from a potential) , and it leads to values of J„ expanded in 
powers of /, together with the determination of the successive c„. 



According to the rotation law adopted, the following results are obtained. 

1. For uniform (rigid) rotation, Chandrasekhar's computations yield a surface 
flatness / that depends on the rotation rate : 



where a = fl^RQ/g. The ratio of central to mean density is Pc/ Pm and g 
is the surface gravity. For the Sun, one adopts the following values : 
n = 461.6 nHz, g = 2.7 W cm/sec^, Rq = 6.96 10^° cm, pm = 
1.409 g/cm^, pc = 90 g/cm^. Accordingly, we flnd a = 2.17 10~^ and / = 
1.11 10-^ 

2. For differential rotation : 

(a) The flatness is decreased from 1.11 10^'' to about 8.85 10~^, in con- 
tradiction with the theoretical work of Maeder (1999). This can be 



4 Results 



/ = (0.5 + 0.856pm/pc)a 



(11) 
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explained only if dfl/dr is > 0, a behaviour observed at latitudes 
greater than 50° (Corbard, 2000; Basu et al, 1999). It is not yet 
clearly established whether this flatness is time-dependent or not. 

(b) Different values of the successive gravitational moments, theoretical 
or deduced from observations, are given in Table 

It can be seen that the first method, using helioseimic data, leads to multipole 
moment values lower than those obtained with the second method. However, the 
octopole moment, J4, is much more sensitive than the quadrupole moment, J2, 
to the presence of latitudinal and radial rotations in the convective zone. Taking 
into account recent values of the subsurface rotational gradient, as reported 
by Corbard and Thompson (2002), the two moments J2 and J4 are extremely 
sensitive to physical mechanisms acting just below the surface. This point seems 
to confirm the existence of a very thin transition layer, as we already suspected, 
that we call the leptocline (Godier and Rozelot, 2001). 



5 Relevance of precise solar gravitational moments 

The determination of the outer shape of the Sun has at least four major 
astrophysical consequences : 

1. Accurate determination of the successive gravitational multipole moments, 
and J2 in particular. Until recently, J2 was not really known, values ranging 
from 1.0 10"^ (Ulrich and Hawkings, 1981) to 36 10"'^ (Gough, 1982). How- 
ever, on one hand, it is expected that the fiattening / be lower than Roche's 
limit, which is, in the solar case 2.7 10~^. This means that J2 must be lower 
than 1.18 10~^. On another hand, we were able to fix an upper bound to the 
quadrupolar moment, J2 max = 3.0 10~^, meaning / is lower than 1.38 10~^ 
(Rozelot and Bois, 1998). Beyond this value, there would be observational in- 
compatibilities for specific astrophysical phenomena such as lunar librations. 
Helioseismic data always lead to lower values of the calculated parameters. That 
is the case for instance for the photospheric radius (6.95508 10^" cm instead of 
6.9599 10^° cm), for the fiattening itself (6 10"^ instead of 9 10~^), for the veloc- 
ity rate at the equator (452.3 nHz instead of 467.3 nHz (Javaraiah, 1995^, and 
also for J2. In this context, Pireaux and Rozelot (2003) adopted the following 
theoretical range for the solar quadrupole moment, based on helioseismic data 
(Paterno et al., 1996 ; Pijpers, 1998) and on a solar stratified model taking into 
account latitudinal differential rotation (Godier and Rozelot, 1999), 

J2 = (2.0 ± 0.4) lO"'^ (12) 
whereas the theory of Figures (Lefebvre, 2003) yields 

J2 = (6.52 ±2.55) 10"^ (13) 

■'solar rotation rate deduced from the sunspots is greater than 461 nHz, see for instance 
Table 3.2 in Lefebvre, 2003. 



9 



These discrepancies still need to be resolved. 

2. Accurate determination of the perihelion advance of a planet, for which a 
knowledge of the solar quadrupole moment would help to determine the rel- 
ativistic post-Newtonian parameters 7 and p. Indeed, the perihelion advance 
(with respect to the classical Keplerian prediction) of planets (AW) is a combi- 
nation of a purely relativistic effect and a contribution from the Solar quadrupole 
moment. It is given by the following expression"' : 



AW = AWo.gr 



\{2a^ + 2aj~(3)-F(a)J2 



(14) 



where AWq^gr represents the contribution due to General Relativity with J2 set 
to zero ; a, /?, 7, the parameters describing the theory of the gravitation in the 
parameterized post-Newtonian formalism ; T{a) is a coefficient depending on 
the considered planet (see Pireaux and Rozelot, 2003) for further developments. 
In General Relativity, Q! = /3 = 7=:1, so that the perihelion advance of planets 
can be written as : AW = 42.981 [ 1 - J" J2] (arcsec/cy), with .F = [RQC^iSsin^i- 
l)]/[2G'Moa(l — e^)] where G, c, Rq and Mq have their usual meaning ; and a, 
e and i are respectively the semi-major axis, the eccentricity and the inclination 
of the orbit of the body in question. The contribution of J2 to AW is relevant. 
In the case of Mercury, J^= +2.8218 10^. 

The parameter 7 is constrained by the CASSINI doppler experiment (Bertotti 
et al, 2003) to 

7- 1 = (-2.1 ±2.3) 10"^ (15) 

while observational constraints on the Nordtvedt effect (Williams et al., 2001) , 
combined with the observational range for 7, leads to 

(3 e [0.9998; 1.0003] (16) 

Hence, it turns out that General Relativity is not excluded by those Solar Sys- 
tem experiments. However, General Relativity would be incompatible with the 
Mercury perihelion advance test if J2 = was assumed. But with the adopted 
theoretical range for J2 in Eq. General Relativity agrees with this latter 
test... and there is still room for an alternative theory too (see Fig.|2||. 

3. Accurate determination of the precession of the orbital plane of a planet 
about the Sun's polar axis. J2 also contributes to this relativistic effect, which is 
unaffected by the relativistic PN combination {2a^ + 2aj — (3). The precession 
of the orbital plane is more easily discernible for moderately large values of the 
inclination i. 

4. Modern planetary ephemerides now include a non null value of J2. However, 
they are presently not able to infer simultaneously all the usual parameters 
(masses, radii,..., PN parameters) and J2 from fits to observational data, due 
to strong correlations. The order of magnitude adopted for J2 corresponds to 

■^In some references, the coefficient of the term containing the contribution of the orbital 
inclination is improperly written. 
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Fig. 3 — For a given value of J2 the perihelion advance of Mercury constitutes a test 
of the PN parameters /3 and 7. 

In the (3 fEa. ll6|l and 7 fEa. ll5|l plane (a set to 1), we have plotted la (the smallest), 
2(7 and 3o" (the largest) confidence level ellipses. Those are based on the values for the 
observed perihelion advance of Mercury, AWobs, given in the literature and summa- 
rized in reference (Pireaux and Rozelot, 2003). Remark also that the position of the 
ellipses varies according to the value of J2 chosen ; but, their orientation is determined 
by the combination {2a^ + 2aj — (3) that appears in the expression for AW . Never- 
theless, General Relativity is still in the 3(t contours for the allowed theoretical values 
of J2 argued by the authors fEn. ll2ll . Fig. 1 a, b, c represent the confidence contours 
for /3, and 7, J2 fixed to its minimum, average and maximum value respectively (Eq. 

10"'', but, for example, the estimation of the PN parameters is rather tolerant 
to the assumed value of J2. 

6 Observations 

Measurements of J„ are not directly accessible. 
The multipole moments could be measured dynamically by sending and accu- 
rately tracking a probe carrying a drag-free guidance system to within a few 
radii of the solar center. The J„ are then inferred from the precise determina- 
tion of the trajectory (Pireaux and Rozelot, 2003). 

Alternatively, the Jn can be inferred from in orbit measurement of solar prop- 
erties (Rozelot and Godier, 2002). Indeed, solar asphericities, encoded mainly 
by the first two coefRcients C2 and a (Eq. EJ can be observed. An estimate of 
these two coefficients has been recently derived from SOHO-Af-D/ space-based 
observations (Armstrong and Kuhn, 1999) : 

C2 = (-5.27 ± 0.38) 10-6 and C4 = (+1.3 ± 0.51) lO^^. 
These results were obtained by measuring small displacements of the solar-limb 
darkening function (for further details, see Kuhn et al, 1998), and the c„ co- 
efficients are comparable to an isodensity surface level (see footnote 1). >From 
Earth-based observations at the scanning heliometer of the Pic du Midi Obser- 
vatory, we also obtained estimates of C2 and C4 coefficients. We have shown that 
asphericities can be observed. A bulge appears to extend from the equator to the 
edge of the royal zones (about 50° latitude), with a depression beyond (Rozelot 
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Fig. 4 - Observations of the solar asphericities by means of the scanning heliometer 
at the Pic du Midi Observatory, in 2001 (left) and in 2002 (right). On the left plot, C2 
and C4 derive directly from Eq. |2l whereas, on the right plot, these c„ come from a fit 
to the best ellipsoid passing through the measured points. 



et al, 2003). Such a distorted shape can be interpreted through the combination 
of the quadrupole and hexadecapole terms, which, as shown previously, directly 
reflect the non uniform velocity rate in surface. Moreover, this distribution im- 
plies a thermal wind effect, blowing in the low density surface layers of the Sun, 
from the poles towards the equator, just as in terrestrial meteorology. This phe- 
nomenon was recently studied by Lefebvre and Rozelot (2004) who explained 
for the first time the existence of small distortions, not exceeding some 20 mas 
of amplitude. Fig. ^ illustrates those observations. Accordingly, We found : 

Year 2000 : C2 = (-7.6 ± 0.2) 10"^ and C4 = +2.2 lO''^ 
Year 2001 : C2 = (-1.1 ± 0.5) 10"^ and C4 = +3.4 10"^ 
Year 2002 : = (-3.8 ± 0.8) 10^^ and C4 = +2.5 10"^ 

The differences in the estimates stem from the difficulty of observations, mainly 
due to seeing conditions. The mean observed value of C2, -7.5 10~^, is not too 
far from the theoretical one, « -(2/3)/«-5.9 10~^ with / = 8.9 10"^ based on a 
solar model with a differential rotation law. The coefficient C4 remains difficult 
to match with the theory, which predicts +(12/35)/^ for a uniform rotation law. 
The only explanation is that the distorted shape coefficient C4 is very sensitive 
to surface phenomena. 



7 Conclusions 

The main conclusions that can be drawn from this review are the following : 

1. HeHoseismic data, even when reliable, imply a solar surface rotation rate in 
conflict with other observed photospheric rotation data. As a consequence, the 
values of both the solar gravitational moments ( J„) and the shape coefficients 
(c„) are slightly different, larger in the latter case. These discrepancies still need 
to be resolved. 

2. The latitudinal dependence conveys sub-surface physical mechanisms that 



12 



can be explained theoretically. Thus, in spite of the fact that ground-based 
observations are altered by seeing effects amplifying or superimposing noise, 
it could be suggested that the solar shape is not a pure spheroid. It is clear 
that improvements on this question will depend on better understanding of the 
sub-surface layers and on better measurements of the limb shape. 

3. The exact shape of the Sun critically depends on the rotation law in the 
external layers. 

4. In the setting of General Relativity, it is clear that accurate measurements 
of the perihelion advance of small planets such as Icarus will help to better 
determine the quadrupole moment (dynamical estimate). Measurements of the 
multipole moments ( J„) will be obtained as a by-product of two dictinct space 
missions : Beppi-Colombo (scheduled for 2009) and GAIA (scheduled for 2010). 
These missions are essential for future developments on that question. 

0. In-orbit solar measurements with sufficient resolution are also essential. 
Measurements of the asphericities (c„) and determination of the true helioid 
are one of the main goals of the future space mission PICARD, scheduled to 
be launched by 2008 ; and also the main targets of balloon flights, such as SDS 
flights. 
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Note 1 : SDS stands for Solar Disk Sextant 

Note 2 : SSE stands for Stellar Structure Equations (based on Ens.l4land l6^ . 

Note 3 : The apparent large error comes from the fact that the value is a weighted average of several rotation rates. 

Note 4 : A mistake has been made in the second term of J4 (i,e. mA^), which was incorrectly multiplied 



by "/" in the computations. 

Tab. 1 — Solar gravitational multipole moments quoted from different authors and 
methods. 
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